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In recent years, a number of nodal superconductors have been identified; d-wave superconductors 
in high T c cuprates, CeCoIns, and k-(ET)2Cu(NCS)2, 2D /-wave superconductor in Sr 2 Ru04 and 
hybrid s+g-wave superconductor in YNi2B2C. In this work we conduct a theoretical study of nodal 
superconductors in the presence of supercurrent. For simplicity, we limit ourselves to d-wave and 
2D /-wave superconductors. We compute the quasiparticle density of states and the temperature 
dependence of the depairing critical current in nodal superconductors, both of which are accessible 
experimentally. 

PACS numbers: 74.25.Sv, 74.20.Rp, 74.25.Fy 



I. INTRODUCTION 

Since the discovery of the organic (TMTSF) 2 PF 6 and 
heavy fermion CeCu2Si2 superconductors in 1979, the 
high T c cuprate superconductors LaBaCu04 in 1986 and 
Sr 2 Ru04 and YNi 2 B 2 C in 1994, the new class of nodal 
superconductors took center stagei However, until re- 
cently, except for <i-wave symmetry in high T c cuprate 
superconductors^^ the gap symmetry has not been ex- 
plored directly for other nodal superconductors. 

In the last few years, it has been recognized that the 
Dopplcr shift in the quasiparticle spectrum in the vortex 
state 4,5 provides extremely sensitive means of studying 
the nodal lines and points in the energy gap A(k)MiZ&& 
More recently the above analysis has been extended 
to the hybrid s+g-wave superconductors— Indeed, 
through angular dependent magnetothermal conductiv- 
ity measurements on high quality single crystals at low 
temperatures (T < IK), Izawa et al. have identified 2D 
/-wave superconductivity in Sr 2 Ru04^ (i-wave super- 
conductivity in CeCoIn 5 14 and K-(ET) 2 Cu(NCS) 2f ±£ as 
well as s+g-wave superconductivity in YNi 2 B 2 C— 

However, as seen in Ref.llTt the nature of superconduc- 
tivity in Sr 2 RuC>4 is still controversial. Clearly both the 
specific heat data and the magnetic penetration depth 
data indicate the presence of line nodes in Sr 2 RuC>4. 8 ' 18 ' 19 
Further, both the angular dependent thermal conductiv- 
ity data— and the ultrasonic attenuation data 20 exclude 
models with vertical line nodes— *2£ Finally the ther- 
mal conductivity data by Izawa et ali^ and Suzuki et 
al£L are incompatible with the two gap model by Zhit- 
omirsky and Rice— — Therefore it is important to un- 
derstand whether /-wave superconductivity is realized in 
Sr 2 RuC>4. More recent angular dependent magnetospe- 
cific heat data for T > 0.1 5K by Deguchi et aZ— do 
not change this situation, though some of the horizontal 
nodes may be converted into a set of point minigaps. 27 

In this paper, we concentrate on two superconducting 
order parameters^ 

A(k) = Acos(20) d-wave, (1) 

A(k) = dAexp(±i(p) cos(x) /-wave, (2) 



where <f> is the in-plane angle measured from the a-axis 
((f) = tan -1 (ky/k x )), and x — c kz with c, being the o 
axis lattice constant. We study the quasiparticle density 
of states (DOS) in the presence of superconductivity (SC) 
and external current, which is accessible through tunnel- 
ing spectroscopy. We also examine the dependence of 
critical current in the superconductor as a function of 
temperature. 

Making use of the mean field quasiparticle Green func- 
tion we can express the quasiparticle density of states 
in the presence of uniform supercurrent in terms of 
simple integrals which are evaluated numerically. The 
superconducting order parameter is also determined 
within the mean field theory. Then, as in the s-wave 
superconductor^ the order parameter A(T) is modified 
in the presence of supercurrent. We see that the super- 
current first increases linearly with the pair momentum 
q s , reaches a maximum value, and drops as q s is further 
increased. We call this maximum value the depairing 
critical current in contrast to the usual critical current 
associated with depinning of vortices and vortex lattices. 
The advantage of studying the supercurrent, compared 
to other transport properties, lies in its sensitivity to the 
condensate itself. 



II. QUASIPARTICLE DENSITY OF STATES 

In the presence of superflow, the quasiparticle Green 
function in the Nambu-Gor'kov formalism is given 

Cr -1 (iWn, k) = iuj n - v F ■ q s + £ k p 3 + A(k)pi<7i, (3) 

where <7i and pi are Pauli matrices acting on the spin and 
Nambu indices respectively. In the quasi-2D system we 
are considering 31 

ik = ^-{kl + k 2 y )-2tcos X - P (4) 
and q s is the pair momentum. 
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The quasiparticle density of states in the presence of 
supercurrent is given by 



f-wave, J || a or d-wave, J || c 



iV(£) = -±VlmTr[G(£,k)], 

k 

which can be reduced to 
g(E) = N{E)/N = Rei 



(5) 



\E-v F - q s 



y/(E-v F •q fl )2-AV 2 / 

(6) 

where iVo = 2m /tt, the function / is cos(2r/>) for the <i- 
wave and cos(x) for the /-wave cases, and the averaging 
is done over angles (•••)= ( 2 n) 2 f d<j> f dx (■ ■ •) ■ For sim- 
plicity we shall consider the following four cases. 



A. /-wave SC with J || a and d-wave SC with J \\ c 

The quasiparticle density of states (DOS) is given by 

(7) 



i r 

g{E,s) = - / d<j>g d (E,s,(p), 



where s = vf<1s for in-plane current or s — 2tq s = vfcQs 
for c-axis current, and 



2 f / A 

s, 0) = - Re \K 



\E — scosi 



(8) 



The quasiparticle density of states in the absence of su- 
percurrent for both /-wave and d-wave cases is given by 
g d (E, 0, 0)M Here K(k) is the Complete Elliptic Integral 
of the First Kind as a function of the modulus^ k and 
Re{K(k)} = k^Klk- 1 ) for k > 1. 

The quasiparticle DOS is shown in Fig. ^ I n the ab- 
sence of current, say for the ci-wave case, the quasiparticle 
spectrum has a line of degenerate saddle points along x 
which give rise to a logarithmic peak in the DOS. At 
finite current, the Doppler shift s smx breaks this degen- 
eracy, leaving only two discrete saddle points. Thus the 
logarithmic peak at |J5| = A is split into two cusps at 
|£|=A±s. 



B. /-wave SC with J || c 

The integral reduces to 



\E - ssinxl 



1 /-t/2 

g(E, s) = - Re / d x ^=- : , _. 

77 J-tt/2 y [E — s sin X) cos x 



(9) 

In the limit E — > 0, g(0, s) = ^ s2 s +A2 ■ The quasiparticle 
DOS is shown in Fig. The Doppler shift ssinx does 
not lift the degeneracy of the line of saddle points in 
the quasiparticle spectrum along <fi, but does shift them. 
Thus the logarithmic peak at \E\ = A is shifted to \E\ = 
VA 2 + s 2 . Also, the DOS flattens out for \E\ < s as 
more low energy states are added near the gap nodes. 
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FIG. 1: Quasiparticle DOS g{E,s) given in Eq. 0. 
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FIG. 2: Quasiparticle DOS g(E,s) given in Eq. Q. 

C. d-wave SC with J || a 

The DOS becomes 



g(E, s) = i Re 



IE 1 — scos< 



tt Jo \/{E — scos0) 2 — A 2 cos 2 (20) 

' (io) 

When £ = 0, we find 



7(0, s) 



s 

7rA 









1 


Rc < 

















(11) 



where 6 = sVs 2 + 8A 2 /(4A 2 ) and a = (1 + b)/2 - 
s 2 /(8A 2 ). 

The quasiparticle DOS is shown in Fig. Just as in 
the previous section, the logarithmic peak at \E\ = A 
is shifted to \E\ = y / A 2 +s 2 /4, while secondary peaks 
appear at \E\ — A ± s due to the interference of the gap 
Acos(2r/>) and Doppler shift scos0 modulations. 
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d-wave, J II a 



2.5 



- 1 5 



1 



0.5 





s/A = 0.0 




,, s/A = 0.2 -- 




s/A = 0.4 




II S/A = 0.6 











0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

E/A 

FIG. 3: Quasiparticle DOS g(E,s) given in Eq. fTUI . 
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FIG. 4: Quasiparticle DOS g(E,s) given in Eq. ifT^ . 



D. d-wave SC with J || [110] 



The DOS is given by 



g(E, s) = - Re 

7T 



IE 1 — s cos < 



(_E - s cos 0) 2 - A 2 sin 2 (20) 

(12) 



For E — 0, we find 



2 f /2A 
5 (0, S ) = -Rejif f — 



(13) 



The quasiparticle DOS is shown in Fig. As in the 
previous section, logarithmic singularities persist, but the 
main peak at \E\ = A is split into two smaller ones ap- 
proximately at \E\ = A ± s/V2 + 0((s/A) 2 ) for small 
s/A. 



III. GAP EQUATION AND CRITICAL 
CURRENT 



Within BCS theory, the gap equation is given by 

A(k) =T^^y kp Tr[ /9l <7 1 G(ia; n ,p)] ) (14) 



where A(k) = A/(k) and V kp 
tion can be reduced to 



Vf(k)f(p). Thisequa- 



1 



hi 



A_ 
A^ 



= Re ( f 2 cosh 



sz \ 

a/; 



(15) 



at T = 0, where z = cos(0 — (f>o) or z — sinx for the 
in-plane or c-axis current respectively, with 4>o being the 
direction of the in-plane current. Ao is the order param- 
eter at T — and q s = 0. At zero temperature, the 
/-wave case was previously investigated in Ref. [2^, while 
the s and c?-wave cases were considered in Ref. 41. 

The dependence of the order parameter on s is shown 
in Fig. [S] Unlike in the 3D s-wave superconductor, both 
the order parameter A and the supercurrent j s jump (see 
Fig. El at s/Ao ~ 0.8, except for the case of an /-wave 
superconductor with c-axis current. The values of A ob- 
tained via the gap equation and j s are shown for s in- 
creased from zero. In all cases the jump disappears at a 
finite temperature and the transition to normal state be- 
comes continuous. These temperatures are 0.06A and 
0.25Ao for the d-wave superconductor with current along 
the a and [110] axes respectively. The transition is con- 
tinuous for all T > in the /-wave case with in-plane 
current. 

This behavior is strongly reminiscent of the case 
of Pauli paramagnetism in both the s-wave and d- 
wave superconductors^*^ These jumps induced by the 
Pauli term signal the presence of more stable inho- 
mogeneous superconductivity, the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) statei*^*^ Very recently, pres- 
ence of the FFLO state has been reported from 
thermodynamic, 36 ultrasonic^ and temperature depen- 
dent upper critical held measurements^ in ci-wave super- 
conducting CeCoIn5>i4 

Therefore, we speculate the appearance of inhomoge- 
neous superconductivity in the vicinity of s/Aq ~ 0.8 
similar to the FFLO state, but generated by uniform su- 
percurrent, in all the cases we have considered except for 
an /-wave superconductor with c-axis current. A more 
detailed analysis of the free energy and possible FFLO 
state will be reported in the near future. 

The finite temperature gap equation reduces to 




: Re (/ 2 cosh" 1 (| 



2 r ^h(E + sz) + h(E- S z)\ 

in / ' v ' 



A|/| 



VE 2 - A 2 / 5 



where h(e) = np(e) — uf(£, T = 0), with np(e) being the 
Fermi distribution. 
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FIG. 5: Zero temperature order parameter as function of s, 
as defined by Eq. 1151 . 



FIG. 6: Zero temperature supercurrent as a function of s, 
given in Eq. (1181 . Note that the Jo differs by a factor of 
vfc/vf between the in-plane and c-axis cases. 



Finally, we study the behavior of the supercurrent, 
which can be evaluated as follows 



neq s 



and reduces to 



•T^^ev k Tr[G(w n) k)], (17) 



J s 1 -2Re 




A 2 / 2 
(sz) 2 



(18) 



at T = 0, where Jo is a constant differing by a factor 
of vf c /vf between the in-plane and c-axis cases. Solv- 
ing Eq. llrJ|) for given s, and substituting the result into 
Eq. H18|) . we can also evaluate j s for given s. The be- 
havior of the supercurrent as a function of s is shown 
in Fig. |S| For the /-wave case, an analytic expression 
for the zero temperature supercurrent and the depairing 
current /max (at which dj s /ds = 0) have been evaluated 
in Ref. |2JJ. Also, the current curves for the s and d- 
wave cases were numerically evaluated in Ref. |4l[ As 
expected, for large s the equilibrium value of the current 
is zero since superconductivity is destroyed. For small 
s, the current rises linearly, it then reaches a maximum, 
and finally decreases to zero in the same way as the order 
parameter. 

The finite temperature expression for the current be- 



Jos 1 



A 2 f 2 

2Re( z 2 \/l 



+ 2< 



A|/| 



(h(E + sz) - h{E - sz)) E dE 
VE 2 - A 2 / 2 ; 



(19) 



The dependence of the depairing current j max on tem- 
perature is illustrated in Fig. [7\ The behavior of the 
critical current as a function of temperature is qualita- 
tively similar for all of the considered cases except for the 
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FIG. 7: Depairing critical current as a function of tempera- 
ture, defined by the maximum value of j 3 (Eq. 1191 ) for given 
T. 



/-wave SC with current along the c-axis. This case is set 
apart in the same way as when the behavior of the order 
parameter was considered. 



IV. CONCLUSION AND SUMMARY 

We consider the supercurrent as a simple pair breaking 
perturbation which adds a Doppler shift to the quasi- 
particle spectrum. While this study is based on the 
assumption of uniform supercurrent throughout a sam- 
ple, the experimental realization of uniform current is far 
from trivial^ 2 ^ In bulk crystals, the supercurrent will 
be localized near the sample surface due to the Meissner 
effect. Therefore, it is important to have high quality 
thin films or whiskers of relevant compounds. Very re- 
cently, whiskers of Bi2212 have been used to measure 
the cross junction Josephson tunneling)22iiS which im- 
plies that this experiment, while nontrivial, should be 
possible in the future. 
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As long as the condition of uniform current is met, 
the quasiparticle density of states obtained here will be 
readily accessible by scanning tunneling microscopy and 
other tunneling techniques. The temperature depen- 
dence of the critical depairing current may also be mea- 
sured. These measurements will provide useful informa- 
tion on the pairing symmetry of nodal superconductors. 

As we have seen, the DOS is very sensitive to the sym- 
metry of the gap and the direction of supercurrent. In 
principle, the DOS is sensitive to the location and type 
of critical points in the quasiparticle dispersion. These 
critical points are easily perturbed or displaced in differ- 
ent ways by adding to the dispersion an appropriate term 
linear in electron velocity — the Dopplcr shift. Thus, with 
the addition of a uniform current, the change in the sin- 
gularities of the DOS can reveal a lot about the symmetry 
of the unperturbed quasiparticle dispersion. 

Also, in many cases, the spatially homogeneous config- 
uration as considered here becomes unstable for s/A <~ 
0.8. This opens up the possibility of a current induced 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. This 



possibility will be further explored. 

In summary, we have evaluated the effect of the pres- 
ence of finite supercurrent on the quasiparticle DOS of 
nodal superconductors. The quasiparticle DOS can be 
probed by tunneling spectroscopy and is sensitive enough 
to the order parameter and supercurrent direction to 
help differentiate between superconductors of different 
gap symmetries. In addition, we showed the dependence 
of the supercurrent on the external current and the de- 
pendence of the critical depairing current on tempera- 
ture. 
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